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The possibility is investigated that competition between fluctuations at different symmetry-related 
ordering wave vectors may affect the quantum phase transition between a fermi liquid and a lon- 
gitudinal spin density wave state, in particular giving rise to an intermediate 'nematic' state with 
broken rotational symmetry but unbroken translational symmetry. At the marginal dimension the 
nematic transition is found to be preempted by a first order transition but a weak symmetry breaking 
field restores a second order magnetic transition with an intermediate regime in which correlations 
substantially enhance the broken rotational symmetry. Comparison to recent experiments is made. 
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'Stripe' spin densi ty w ave order occurs in many high- 
T c cuprate material s ' 1 * 2 * 3 !. A spin density wave is a lon- 
gitudinal modulation of the spin density 5*(R) char- 
acterized by a wavevector Q giving the periodicity of 
the spin modulation. In a 'stripe', 2Q is not a re- 
ciprocal lattice vector so the magnitude of the spin, 
as well as its direction, varies from lattice site to lat- 
tice site. If the underlying lattice has sufficient sym- 
metry, stripe ordering may occur at one of several in- 
equivalent wavevectors Q a . In the hole-doped high-T c 
cuprates the important physics is two dimensional and 
the lattice has (to a good approximation) square sym- 
metry. For dopings greater than about x = 0.05 and 
less than a material-dependent number ranging from 0.08 
(in YBa 2 Cu 3 6+ ®) to 0.24 (in Nd 2 - x Sr x CuOj^) or- 
der is believed to occur at one of the two wave vectors 
Q x = (tt — S, 7r) or Q y = (7r,7r — 5) wi th 5 doping- 
dependent and generically nonzerc & * 2 * 4 * 7 * 8 * 9 l 

In a stripe state the expectation value of the spin den- 
sity at position R, S(R), may be written 

(£(R))=i a cos(Q a -R + a ) (1) 

The state defined by Eq [l] breaks spin rotation, lattice 
translation and lattice rotation symmetries. A phase 
characterized by a nonvanishing < 5(R) > is also char- 
acterized by nonvanishing 

(T(R)) = (s(R) • 5(R)) - ((s(R) ■ 5(R))) 

~ cos(2Q a • R + 26» a ) (2) 

where the double-bracket indicates also an average over 
position. T is spin rotation invariant but breaks lattice 
translation and rotation symmetry if 2Q is not a recip- 
rocal lattice vector. (One may also consider bond order 
involving < 5(R) • S(R' ^ R) > but this will not be 
important here). T couples linearly to lattice distortions 
and the electronic c harge d ensity, so is observable in scat- 
tering measurements^^ and is sometimes referred to as 
'charge order'. 



The 'stripe' state is also characterized by nonvanishing 

(r?(R)) = ((s Qx (R) ■ S Qx (R) - Sq v (R) ■ S Qu (R))) 

(3) 

where SQ a indicates spin fluctuations with wave vectors 
near Q a . < r\ > is invariant under spin rotations and 
lattice translations but breaks the discrete lattice rota- 
tion sym metry and may be referred to as a nematic order 
parameteil^iSI 

The three bro ken symmetries may be restored at sep- 
arate transitions^^. (Very similar phenomena are well 
understood in the classical physics context of smetic and 
nematic liquid crystals^). If spin order is destroyed by 
fluctuations in the direction of A (as would happen in 
a two dimensional model with Heisenberg symmetry at 
any T > 0), < T > and < rj > may be expected to re- 
main non-zero. If long ranged order in T is destroyed by 
fluctuations in 9, < rj > may remain nonvanishing. In 
physical terms, the state with < S >=< T >= but 
< r\ >=/= has the property that fluctuations around one 
of the Q a are larger than fluctuations around the other 
possible Qb^a- 

Experimental evidence suggests that this sequence of 
transitions indeed o ccur s in some high T c compounds. 
Many measurements 1 ^ indicate that in underdoped 
cuprates the ground state (if superconductivity is sup- 
pressed) is characterized by magnetic scattering at the 
wavevectors Q x ,y and 2Q XiJ/ but not at Q x ±Q y , imply- 
ing that the scattering signal arises from a superposition 
of domains with order at either or Q y . As tempera- 
ture is raised above an ordered state the Bragg scattering 
at Q vanishes first, leaving an intermediate state with 
Bragg scattering only at 2Q X and 2Q y H Recent neu- 
tron scattering measurements 4 ^ on a monodomain sample 
of Y Ba2Cu30eA5 indicate a wide temperature regime 
where there is order neither at Q x-y nor at 2Q x ^ y but 
where the fluctuations associated with ordering wavevec- 
tor Q x have much longer spatial range and stronger 
temperature dependence than the fluctuations associ- 
ated with potential ordering wave vector Q y . Trans- 
port measurements have detected rotational symmetry 
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breaking and, recently, an enhancement of the N ernst 
effect in this temperature regime has been report e d 9 * 12 *, 
also consistent!^ with an intermediate nematic phase. 
Similar transport behavior in the (Nd/ Sr)2CuO^ fam- 
ily of materials has also been interpreted in terms of a 
nematic phase or regime 6 -. It is however important to 
note that the crystal structure of both of these materi- 
als is such that a CuOi plane is orthorhombically dis- 
torted so it may be more appropriate to describe the 
observed 'nematic' regime as being characterized by a 
strong and strongly doping and temperature-dependent 
enhancement of a pre-existing anisotropy. 

Closely related issues hav e been discussed in the con- 
text of the pnictide material J 14 * 15 ! where 2Q is a recipro- 
cal lattice vectoi!^ so only the spin and nematic orders 
arc relevant. Also in pnictides a strong coupling to lattice 
distortions believed to be important. 

These and related experiments have focussed theoreti- 
cal attention on 'nematicity'. A number of works consider 
nematic phases which are taken to b e conceptually inde- 
pendent of any density wave orderinj- 13 | 17 | 18 | 19 | 20 | 21 | 22 | 23 l 
This paper considers the density wave instability as pri- 
mary with the nematic phase arising from it. The phys- 
ical idea is straightforward: in a stripe situation, den- 
sity wave ordering at one possible ordering wavevector 
Q Q must act to suppress density wave ordering at the 
other possible wavevectors Q&^ a - Thus, as a putative 
'stripe' quantum critical point is approached competition 
between fluctuations at different wavevectors may drive 
a 'nematic' transition at which the system chooses one 
wavevector at which the fluctuations will become criti- 
cal, while fluctuations at the other wavevectors remain 
massive. Alternative possibilities are that the compe- 
tition is important only for selecting the relevant state 
inside the density wave ordered phase, or that competi- 
tion between fluctuations may drive the transition first 
order. In renormalization group language the question is 
whether there is a relevant operator at the stripe critical 
point and, if so, does it imply a flow to a new 'nematic' 
critical point or a runaway flow indicating a first order 
transition. 

This paper approaches the physics in terms of a T = 
instability of a disordered fermi-liquid phase using the 
standard 'Hertz' model of a densit y wave transition in 
a two-dimensional fermi liquicP^H. The main finding is 
the phase diagram depicted in Fig [TJ for lattices with 
square symmetry the quantum 'nematic' transition is 
typically preempted by a strongly first-order transition 
directly to a density wave state; however, a weak explicit 
symmetry breaking restores a continuous transition. 

Incommensurate density wave t ransitions are the sub- 
ject of an extensive literatur e! 26 * 27 !, but the issue of inter- 
est here has been less studied. Physics similar to that 
of interest here has been explored in the context of clas- 
sical spin models for cuprates^ and pnictides^. Also, 
although their main focus was on transitions between 
nematic and density wave ordered states, Sun and co- 
workers observed^ that the basic fermi liquid to den- 
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FIG. 1: Zero temperature phase diagram of two dimensional 
stripe ordering model in plane of doping 8 and tetragonal 
symmetry breaking parameter S. Solid lines: phase transition 
between paramagnetic metal ('disordered') and density wave 
("stripe") state. Star (light blue on-line): tricritical point sep- 
arating first order (heavy line) and second order (light line) 
phase transitions. Dashed line (black on-line); putative sec- 
ond order density wave transition preempted by first order 
transition. Small filled circles (red on-line): paramagnetic ne- 
matic phase, defined only at S — and separated from density 
wave and disordered phases by second order phase transition 
points marked respectively by large circle (green on-line) and 
square (purple on-line). Nematic phase and associated tran- 
sitions are also preempted by first order transition 



sity wave transition would likely be first order. DePrato, 
Pelisseto and VicarP^I used renormalization group tech- 
niques to classify the quantum critical fixed points of a 
model of 'stripe' quantum criticality involving undamped 
spin excitations in two spatial dimensions, in particular 
identifying and analysing regimes of stable second order 
transitions. At these transitions there would be no inter- 
mediate nematic phase separating the density wave and 
disordered state. 

Pelissetto, Sachdev and VicarPS studied a spin den- 
sity wave transition occurring inside a d-wave supercon- 
ducting state. The most relevant perturbation involved 
coupling of nodal fermions to a nematic order parameter 
derived from the spin fluctuations in a manner very simi- 
lar to what is considered here; however again the density 
wave transition was not preempted by a nematic one. Qi 
and XrP^ studied a spin-fermion model, finding a run- 
away flow also indicating first order transitions, but only 
at an exponentially long length scale. 

A theory of the disordered spin density wave may be 
obtained from Eq [7] by regarding A and 6 as slowly fluc- 
tuating quantities. We specialize to the two dimensional 
square (or rectangular) lattice. Because we are interested 
in the transition from a fermi liquid where all spin ampli- 
tudes are small we combine A a and 9 a into a new complex 
field tf a (r) = A a {r)e i6 "^ which we assume is described 
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by the action S = S, 



dyn 



with 



a— x,y 

d 2 rdr — X— ( ■ i> x 



1p* ■ lp a 



2 l f 



(4) 



v 
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The 5 a are control parameters (for example, doping) 
which tune the system through the magnetic quantum 
critical point and we have allowed for the possibility that 
deviations from tetragonal symmetry favor ordering in 
one direction rather than another. Because the experi- 
mental evidence in high-T c materials indicates that spin 
fluctuations are stronlgy peaked near discrete momen- 
tum values we do not need to consider the possibility of 
a continuous rotation of the wavevector and any potential 
nematic phase would have a strong Ising anisotropy. 

Of the six possible quartic nonlinearities (see^l for a 
complete list) Eq [4] includes only the two which are im- 
portant for the present purpose, neglecting terms which 
favor spiral and other non-stripe states or renormalize 
the basic stiffness against large amplitudes of the fields. 
The crucial term is the third one, which quantifies the 
extent to which orders in the x and y directions compete 
with each other. The relevant case is v > u, so that fluc- 
tuations compete. On the mean field level stability of 
this theory requires that u > and v > —u. 

We assume standard overdamped dynamics. For ease 
of writing we present Sd yn in frequency space: 



S, 



dyn 



' e^e 



a—x,y 



d 2 r^|^-^| (5) 



The theory requires an ultraviolet cutoff. We measure 
energy in units of T and impose a hard cutoff, elimi- 
nating all processes for which |f2|/r + k 2 > A. Qi and 
Xil^ studied essentially this model, but with an addi- 
tional (tpi ■ V>2) 2 coupling. 

As defined the upper critical dimension of the model 
is d = 2 and the physics may be studied by a 
renormalization-group analysis. The required beta func- 
tions are given in Eq (3.1) oP^. It is useful define new 
variables g and <p hy u — gcos(f>, v = gsiacj) which flow 
according to (here the dot denotes changes with renor- 
malization group cutoff parameter) 

g = - (7 cos 3 <j)+ 11 cos sin 2 <j>+ 2sin 3 0) g 2 (6) 
= (3 sin 3 — 2 sin 2 <f> cos — sin cos 2 0) g (7) 

Because g ~ g 2 while ~ g, flows much more raipdly 
than g and the content of the theory may be understood 
from a constant g. In Eq|7]the angle 0/4 (corresponding 
to u = v) is a separatrix. For (j> < ir/4 the flow is to- 
wards = 0, but in the <p > 7r /4 case relevant to stripe 
physics the flow is towards a fixed point value which is 



close to 7r. Noting that u turns negative at = 7r/2 
we see that the renormalization group analysis indicates 
that when the flow passes this point the two dimensional 
stripe fixed point becomes unstable towards a first order 
transition. The basic conclusion is perhaps not surpris- 
ing: the model of two coupled order parameter fields is a 
textbook exa mple of a runaway flow leading to a first or- 
der transitio n 1 32 * 33 ! . De Prato et aPSl observed that their 
more general renormalization group equations had only 
unstable fixed points near the marginal dimension and 
Sun et. al. noted that multicritical points of this type 
tend to be unapproachable due to the presence of run- 
away flows. 

Qi and Xu^ similarly noted the possibility of a run- 
away flow. The equations of^ involve three couplings 
and are thus more complicated to solve. A numerical 
solution was presented which indicated a runaway flow, 
albeit beginning at an exponentially low scale, whereas 
what is found here is a first order transition at a scale 
which is not, in general, exponentially small. 

The first order nature of the transition arises from com- 
petition between fluctuations associated with the two or- 
dering wave vectors. An explicit symmetry breaking term 
(arising e.g. from the chains in YBCO) would grow 
under renormalization and if it became large enough, 
would quench the fluctuations at one of the two wave 
vectors, thereby permitting a continuous behavior. To 
understand the energy scales involved in this scenario 
we consider a self-consistent one-loop analysis, which 
while less rigorous than a renormalization group treat- 
ment has a transparent physical interpretation and allows 
for straightforward estimations of energy scales. 

To implement the self consistent one loop theory we 
write the model as a functional integral and decouple the 
nonlinearities |i/ , :c | 2 + |V'y| 2 an d \ipx\ 2 — IV'J 2 by Hubbard- 
Stratonovich fields iX and r\ respectively and then inte- 
grate over the ip fields obtaining the action 

S[\,V\ = -r^— + ^^- + lTrHTl + 5 x + i\ + r ] ] 
2{v + u) 2(v — u) 2 



-Trln[n + <5^A-7;] 



(8) 



with r±o = |f2| + k 2 . Eq [8] is written for the param- 
agnetic phase; the factor of 3 is the spin degeneracy. 
Mean field theory corresponds to finding the A* and rf 
which extremize S. The extremal values of A are imag- 
inary; we write iX = 5 — r with 5 = (S x + S y )/2 and 
introduce A = (6 X ~ 8 y )/2 which parametrizes any ex- 
plicit breaking of tetragonal symmetry. The correlation 
length S >x _ y for spin fluctuations around the wavevectors 
Qx, v is CaT. 2 = r ± (77 + A). In the paramagnetic phase 
> 0. If one of the fields, say ip y , orders, then 
the Heisenberg symmetry ensures that the two trans- 
verse components are gapless (£j 2 y = 0) while the longi- 
tudinal component has a correlation length determined 
by the magnetization m, (£^ 2 = rn 2 ) so that the term 

3 
2 



Tr In [I±o 



iX - i]} -> Tr ln[n ] + \Tr ln[II - 
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Fang et. al!^ presented a large-N analysis of a classical 
spin model which leads to a theory very similar to that 
defined by Eq[8]if S is chosen to be deep in the ordered 
phase so quantum fluctuations are unimportant. 

Defining 6 cr it to be the value at which S is extremized 
at A = r x — r y = 0, redefining 5 XtV as the difference 
from 5 cr it, introducing u(v) = 3u(v)/4ir 2 and explicitly 
carrying out the minimization in the paramagnetic phase 
at T — we obtain: 



5x = f x ( 1 + u hi — ) + vr y In — 



(9) 



/•„ I 1 a In —| + vr x ln — (10) 



In the tetragonal symmetry case A = and at 6 > 
Eqs 9|10 admit an isotropic solution r x = r y > corre- 
sponding to the conventional paramagnetic phase. How- 
ever, for v > u one finds that as 8 is decreased below 
a critical value 5 nem > the isotropic solution under- 
goes a bifurcation to a solution with r x ^ r y . (Fang 
et. al.^ found a very similar transition, in their case 
thermally driven). This is a nematic phase: the nematic 
order parameter is 77 ~< \ip x \ 2 > — < \ip y \ 2 >■ Writing 



V 2 , r y 



rj/2 and linearizing in 77 we find 



2v + v 2 



(11) 



However, within mean field theory this critical point is 
typically preempted by a first order transition to a state 
with long-ranged stripe order. The first order transition 
manifests itself as a failure of numerical routines to find 
a solution to Eqs 9|10 as 6 is decreased below a spin- 
odal value greater than S nem , and may also be seen more 
directly. 

We have computed the energy of a magnetized state 
by generalizing Eq [8] as described above. We find that 
at 6 = 5 nern an ordered state with r = ±77 and m > 
exists and has lower energy than the paramagnetic state, 
provided that v is not too large. Thus at some S > 5 nem 
the system will jump from the isotropic paramagnetic 
phase to a phase with long ranged order. The largest 
6 at which an ordered phase may be sustained may be 
ascertained from the 6 at which there is a solution to Eqs 
9|10 with r y — and r x > 0. This 5 = Sist_ order is 



ve 
v ■ 



(12) 



By comparison of energies we find that if 8ist_ order > 
S nem the nematic transition is preempted by a first order 
transition. As the second equality of Eq [T2| shows, for 
v > v c (u) with v c (u = 0) = e-2 w 0.718.., S nem > S order 
so that a second order transition can exist at large v. 
However, in the large v regime in which the transition 
(within the present theory) is second order, substitution 
into the defining equations shows that the renormalized 




FIG. 2: Dependence on control parameter 8 of inverse corre- 
lation length parameters for x (solid line) and y (dashed line) 
wavevectors for different levels of explicit symmetry breaking 
5 X — S y + A with ansiotropy parameters A as indicated and 
interactions u — 0.3 and v — 0.6. 



mass r is of the order of the cutoff, indicating that the sec- 
ond order transition occurs in a regime beyond the range 
of validity of the critical theory. We interpret this result 
as meaning that the nematicity arising in the large v case 
is an intrinsic phenomenon arising from short length scale 
physics and not directly related to the singular magnetic 
critical fluctuations. 

We are now in a position to consider the effects of an 
explicit breaking of the C4 lattice rotational symmetry, 
such as occurs in Y Ba-2.Cu?,0§-\- x and N ' di- x Sr x CuO^. 
This acts analogously to a magnetic field at a ferromag- 
netic transition and, if large enough, will convert the first 
order transition to a continuous one. The small value 
of (8 v t_ order - 5 nem )/8 nem suggests that the symmetry 
breaking field need not be large. Fig [2] presents results 
obtained by solving Eqs |9|10 for S x = S y + A for a repre- 
sentative choice of parameters u = 0.3 and v = 0.6 (im- 
plying 5ist order = .071 and 5 nem = 0.0328). and several 
levels of anisotropy. One sees that for these parameters 
anisotropy greater than about .01 ~ (disorder — S ne m)/^ 
restores a continuous behavior. 

Extending the results to temperatures T > is com- 
plicated by issues including the difficulty of treating the 
classical transitions below their upper critical dimension, 
the interplay between 'nematic' and 2Q ("charge") or- 
der and the possible first-order nature of the transitions. 
A detailed discussion will be given elsewhere but a few 
remarks are in order here. 

One should distinguish two situations: either the 
model has an intrinsic tendency towards nematic behav- 
ior, not directly related to the stripe criticality (in the 
model defined by Eq [4] this would occur if v is large, 
both absolutely and relative to u) or the first order 
transition preempts the T = nematic phase. Mod- 



els with 'intrinsi c' nematicity have been extensively dis- 
cussed elsewhereP 1 17 1 18 1 19 1 20 1 21 ^ 212 ^ and will not be con- 
sidered here. In the case where the T = nematic phase 
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is preempted by a first order transition there remains the 
possibility of nematic behavior at T > 0. 

Indeed the model defined in Eq [4] trivially exhibits ne- 
matic behavior if parameters are chosen so that the model 
is in the magnetically ordered state at temperature T = 
and if charge order is not important (for example because 
2Q is a reciprocal lattice vector). The ordered state must 
select one of the two wave vectors, say Q x . The Heisen- 
berg symmetry and two dimensionality then implies that 
if temperature is slightly raised above the ordered state, 
fluctuations near have a correlation length of the 
'renormalized classical' form £ ~ exp[p s /T] while fluctu- 
ations near Q y have a relatively short and weakly tem- 
perature dependent correlation length, so the resulting 
state is a 'nematic'. The question then is whether this 
'thermal spin nematic' behavior vanishes via a first or a 
second order transition. If a weak interlayer coupling is 
added to the model one may ask if the nematic behavior 
survives at temperatures higher than the three dimen- 
sional ordering temperature (this issue was addressed in 
the classical model of Ref [2]). If charge order is also 
important, one may ask if the nematic behavior vanishes 
at the charge ordering temperature or at a higher tem- 
perature, and what are the orders of the transitions. 

Clearly a transition which is first order at T — re- 
mains first order as the phase boundary is extended to 
T > 0, at least within some distance of the zero temper- 
ature transition, so an extension of the present results 
to T > would suggest a first order thermal transition, 
at least for dopings near where the T = ordered state 
vanishes, but it is not ruled out that the transition may 
become second order as doping is reduced deeper into the 
ordered phase. 

The self-consistent one-loop approximation employed 
above leads to a reentrant behavior. As T is increased 
from zero the first order transition line extends to the 
larger S side of the critical point and the transition is from 
an isotropic paramagnetic phase at low T to a higher T 
'nematic renormalized classical' phase in which the spin 
fluctuations near one of the two ordering wave vectors 
have a very long and rapidly temperature-dependent cor- 
relation length while those near the other have a corre- 
lation length which is relatively short and weakly tem- 
perature dependent. Whether the 'nematic renormalized 
classical' phase is also characterized by long ranged 2Q 
charge order depends on details. As the temperature 
is further increased the charge order (if any) disappears 
and there is a second transition (typically also first-order 
within this approximation) to a fully isotropic phase. The 
first order transition persists as parameters are tuned to 
move deeper into the insulating phase. While the reen- 
trant behavior seems unphysical, and is likely to be an 
artifact of the approximation, the qualitative result of a 
range of zero temperature parameters above which the 
thermal transition is first order is robust. 

The classical-spin results o£^ provide an interesting 
perspective on this issue. As remarked above, this model 
is in essence the self-consistent one-loop approximation 



to the classical limit of the model defined in Eq [8] for 
5 deep within the T — ordered phase and 2Q a re- 
ciprocal lattice vector. The authors of Ref [H] consid- 
ered a regime in which the physics was controlled by 
the 'renormalized classical' divergence of the correlation 
length and moreover chose parameters corresponding to 
an extremely weak tendency towards nematic ordering. 
RefSl defines a dimensionless quantity K/J2, which is 
essentially (v — u)/(v + u) in the notation of this paper. 
This was taken to have the very small value 7.5 x 10 -3 . 
For this parameter value Ref [3] reported a second order 
'nematic' transition which occurred at very low tempera- 
tures deep in the 'renormalized classical' regime and was 
not preempted by a first order transition. Ref [14] did 
not report results for larger K, but solving their equa- 
tions indicates that as K is increased to a value > 0.075 J2 
the thermally driven transition again becomes first order. 
Now 4-7T J2 is essentially the spin stiffness of the magnet- 
ically ordered state which may be thought to grow with 
distance into the ordered phase so if one imagines that 
K ~ (v — u) is fixed to a value which is not too large, the 
increase in the magnetic spin stiffness distance into the 
ordered phase may eventually drive the model into the 
second order transition regime. 

Of course, the reliability of the self-consistent one-loop 
approximation may be questioned, especially for thermal 
phenomena. The analysis oP^ indicates that this anal- 
ysis is essentially equivalent to a renormalization group 
analysis for models above the upper critical dimension, 
but for models which are below the upper critical dimen- 
sion this is of course not the case. The self-consistent 
one loop approximation has the exponents of the Gaus- 
sian model, and for this reason may overestimate the 
tendency towards nematic order. For example, one may 
define a 'nematic susceptibility' via the correlation func- 
tion < {tpl — i/jy) 2 >. This correlator is closely related 
to the energy correlator that defines the specific heat ex- 
ponent. The gaussian model dramatically overestimates 
the divergence of the specific heat near 2D and 3D clas- 
sical transitions and may well similarly overestimate the 
divergence of the 'nematic susceptibility'. 

Despite these difficulties it is interesting to relate the 
picture presented here to data. The essential point is that 
in models in which 'nematicity' is derived from compe- 
tition between density wave ordering at several wavevec- 
tors, the nematic phase may be pre-empted by a first 
order transition, but if the lattice rotational symmetry 
is explicitly broken a more continuous behavior may be 
restored. There is no direct evidence that stripe order 
terminates at a first order quantum critical point in any 
cuprate, although we note that phase coexistence, a typ- 
ical consequence of first order phase boundaries, is com- 
mon in cuprate materials. The strongest evidence in 
favor of a nematic regime comes fr om th e YBCO and 
{Nd/ Sr)2CuO± families of material d 4 * 9 * 12 * where the lat- 
tice symmetry is explicitly broken. It is very tempting to 
argue that the symmetry breaking pushes the system into 
the continuous transition regime identified above. It may 
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also be worthwhile to reexamine the data to determine if 
a local lattice distortion (favoring one or the other wave 
vector) is present or if a hysteresis has been overlooked. 

In summary, this paper has posed the question of the 
existence of an intermediate nematic phase separating a 
stripe and a fermi liquid phase in terms of fluctuation 
corrections to a putative 'stripe' quantum critical point. 
The physics that can lead to a nematic phase also leads 
naturally to a fluctuation-driven first order transition, 
which near the quantum critical point was found to pre- 
empt the nematic phase. The first order transition is not 
inevitable. If parameters were tuned so that the nematic 
transition occurs 'far' from the putative stripe quantum 
critical point so that the 'nematicity' is an intrinsic effect 
and not driven by critical density wave fluctuations, then 
the considerations of this paper are not relevant. 

It was found that the first order transition could be 
converted to second order by a quite small anisotropy. 
The analysis relied on approximations including the 
self-consistent one-loop theory (which is uncontrolled) 



and the Hertz quantum critical theory (which is sub- 
ject to corrections whose nature remains incompletely 
understood^) but the first order behavior discussed here 
follows from relatively general considerations and seems 
likely to be robust. In systems where nematic behavior 
was found a re-examination of experimental data for sig- 
natures of first-order transitions (for example, hysteresis 
and phase coexistence) may be worthwhile. Extension of 
the results presented here to T > 0, to include coupling 
to the lattice, and to other situations, such as the meta- 
magnetic transition in SV3-RU2O7, would be of interest. 
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